Aging dynamics and the topology of inhomogenous networks 
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We study phase ordering on networks and we establish a relation between the exponent a x of the 
aging part of the integrated autoresponse function \a g and the topology of the underlying structures. 
We show that a x > in full generality on networks which are above the lower critical dimension 
o!l, i-e. where the corresponding statistical model has a phase transition at finite temperature. For 
discrete symmetry models on finite ramified structures with T c — 0, which are at the lower critical 
dimension g(l, we show that a x is expected to vanish. We provide numerical results for the physically 
interesting case of the 2 — d percolation cluster at or above the percolation threshold, i.e. at or above 
g(l, and for other networks, showing that the value of a x changes according to our hypothesis. For 
0(Af) models we find that the same picture holds in the large- N limit and that a x only depends on 
the spectral dimension of the network. 



In recent years, networks and graphs have been suc- 
cessfully applied to very different disciplines 0. This 
extreme flexibility is not surprising, since they are the 
most general way of representing a set of elements con- 
nected pairwise by some kind of relation. In social sci- 
ences, computing, psychology and economy, graphs rep- 
resent communication systems, social relationships, bio- 
logical interactions and statistical models of algorithms. 
In physics, chemistry and biology, they are extensively 
used as models for real complex and inhomogeneous 
structures, such as disordered materials, percolation clus- 
ters, glasses, polymers, biomolecules. Networks repre- 
senting real physical structures are constrainted to be 
embeddable in a finite dimensional space and to have 
bounded coordination number, forming an interesting 
class of structures called "physical graphs" @|- 

Basic physical equilibrium properties of models de- 
fined on physical graphs, among which critical behav- 
iors, depend crucially on the topological features of the 
network. In particular, the relation between large scale 
topology and critical properties is well understood on a 
simple class of network, i.e. regular lattices, where the 
existence of phase transitions, the value of critical expo- 
nents and many other equilibrium quantities depend on 
a topological parameter, the Euclidean dimension d of 
the lattice. As far as systems with continuous symmetry 
are concerned, such as 0(N) models, this general pic- 
ture can be extended to networks. Their topology can 
be completely described by the adjacency matrix A^j, 



dimension" d s of the graph. The spectral dimension uni- 
voquely determines the existence of phase transitions Q 
and controls critical behaviour |J| , much in the same way 
as the Euclidean dimension d does on usual translation 
invariant lattices. In systems with a discrete symmetry, 
although a unique topological indicator, analogous to d s , 
is not known, the fundamental role of topological and 
connectivity properties in determiningequilibrium and 
critical behavior has been pointed out p(. 

The picture is much less clear for out of equilibrium 
processes. Although it is known that dynamical processes 
are influenced by the topology of the network, a direct 
description of this relation is largely incomplete. This is 
the case for the phase ordering dynamics, occurring in 
systems gradually evolving from an homogeneous phase 
to a state with two or more coexisting ordered phases, 
through formation and growth of segregated domains . 
Initally investigated in physics, phase ordering also oc- 
curs in other fields where the topology is of primary im- 
portance, as in biology, where it describes spreading of 
tissues 7], and in economic and social networks, where it 
has been used as a model for propagation of opinions and 
technologies 8] . A deeper understanding of the relation 
between the dynamical process and the topology of the 
network in these frameworks would be of great interest. 

In the simpler case of regular lattices, the influence 
of the Euclidean dimension on dynamical exponents has 
been ascertained in phase ordering in spin systems 0|. 
The generalization of such models to networks is pro- 



with Ai j = 1 if i and j are connected by a link and vided by the Hamiltonian H{a} — J2i j AiJ^i ' °j where 



Aij = otherwise, but the relevant large scale topo- 
logical information is encoded in the Laplacian matrix 
where zi = ^ -Ay is the number of 



A, 



neighbours of site i. Indeed, the density p(X) of eigenval- 
ues of A, p(X) ~ X^'f 2 )- 1 for A -> defines the "spectral 



the indices i and j run over the N sites of the graph, and 
al represents an A/"-dimensional unit vector defined at the 
site i. For Af = 1, H describes the discrete Ising model 
and for Af > 1 the 0(Af) models, with continuous sym- 
metry. Phase ordering dynamics occurs when the system, 



2 



initially at thermal equilibrium at high temperature, is 
quenched to a temperature where the symmetry is broken 
and more phases coexist. In this scenario, an interesting 
physical quantity is the integrated autoresponse function 
x{t, s) describing the effect of a small perturbation acting 
on the system from time s to t: x(t, s) — J R(t,t')dt' , 

where R(t,t') = N~ l £i • ^)/5|^(i')|| Si(t , )=0 is 

the average linear response function associated to an im- 
pulsive magnetic field switched on in i at t' 
hi/\hi\ and (...) means ensemble averages). Interest- 
ingly, x(t, sL can be splitted into a stationary and an 
aging term [<j x(t, s) = Xst{t - s) + x ag {t, s), and the ag- 
ing part features a scaling behaviour with a characteristic 
exponent a x 



Xag(t,s) = S- a xf(t/ S ). 

In phase ordering on lattices [l(J 

f(x) ~ x~ a * for x > 1, 



and a x is related to d by 



d — dj 



(1) 



(2) 



(3) 



for d < du, while a x — 9 for d > djj. Here g?l is the 
lower critical dimension, djj = 3 or djj — 4 for discrete 
and continuos symmetry models respectively, and 9 is the 
exponent regulating the time decay of the topological de- 
fects density p oc t~ e , which does not depend on d. By 
means of Eq. @, a non-equilibrium exponent, a x , is re- 
lated to the topology of the underlying lattice, through 
d. This relation implies that a x > when the system is 
above g?l, i.e. when a phase transition at finite tempera- 
ture T c occurs, a x = at c?l and a x < below dh- 

These results, holding in the simple case of regular 
lattices, suggest that non-equilibrium dynamics of sta- 
tistical models could be related to important topological 
properties of networks. With the aim of testing and in- 
vestigating this hypothesis, we study the response func- 
tion exponent a x on generic physical networks, showing 
that i) a x > is expected in full generality on networks 
supporting phase transitions at finite T c , i.e. above d^', 
ii) for continuous symmetry 0{M) models, in the solu- 
ble large- M limit on graphs the same expression © is 
found, with d s occurring in place of d. This fits with the 
aforementioned general picture that, for continuous sym- 
metry models d s is the topological indicator replacing d 
on graphs; iii) for discrete symmetry models, we provide 
an argument showing that a x = is expected on finitely 
ramified networks (FRNs) 5], which represents the most 
general class of structures without phase transitions at 
finite T c . Although in this case a topological parame- 
ter playing the role of d s is not known and, therefore, 
a generalization of Eq. @ is not straightforward as for 
continuous symmetries, this result conforms with Eq. @ 



and provides, via the existence of phase transitions, a re- 
lation between the positivity of a x and the topological 
feature of the network. Our results are complemented 
by extensive numerical simulations of different discrete 
symmetry models on a representative set of graphs find- 
ing a x — or a x > depending on the absence/presence 
of phase transitions. In particular, we provide results for 
the physically interesting case of the 2 — d percolation 
cluster at or above the percolation threshold, i.e. at or 
above d^, showing that the value of a x changes according 
to our hypothesis. Taking advantage from this, a x could 
be conveniently used to investigate the phase diagram of 
statistical models on general discrete structures. 

Let us start showing that i) a x > is expected for mod- 
els with a phase transition at finite T c . Reparametrizing 
the two time dependence of x(t, s) in terms of the au- 
tocorrelation function C(t,s) = l/JV^ tl (?i(i) • <?i(s)) 
, one obtains the parametric form x(C, s) , which is 
related to the structure of the equilibrium state through 

El 



-T lim 



d 2 x(C,s) 



dC 2 



(4) 



C=q 



where P eq (q)^ = exp[-f3(H{a} + 

H{a'})]6[Q(a,a') — q] is the equilibrium probability 
distribution of the overlaps between two configurations 
{a} and {a'}, Q(a,a') = (l/iV)E^ ' o?, z being the 
partition function. General considerations can be drawn 
on the relation between Xag and the critical behavior 
of the corresponding statistical model. Let us consider 
a system with T c > quenched to < T < T c . In 
this case the form of P eq (q) implies |13| that Xag must 
vanish asymptotically. Then, from Eq. one has 
a x > 0. Even if for quenches to T = Eq. (@J does not 
give informations on Xag, since dynamical exponents are 
known [£j not to depend on T |fj , a x > must be found 
in the whole interval < T < T c . These considerations 
being completely general, they apply to networks as 
well, and a x > for quenches to T < T c in systems with 
finite T c . 

Let us now consider point ii). In the large- Af limit, the 
order parameter </>; can be expanded in the base of the 
eigenvectors of the Laplacian operator and the dynamical 
evolution of the projection 4>k of the order parameter on 
the fc-th eigenvector of the Laplacian operator obeys 0] 



dt 



[A fc + /(t)]0fc + »fc(i) 



(5) 



where is the fc-th eigenvalue of A, I[t) = 

1 / N J2i r ($) - 1 and Vk(t) is the thermal noise. Eq. © 
is formally linear and it can be solved in the large time 
domain. Proceeding as in 0] to compute Xag{t, s) for 
a quench below T c and extracting the leading term, cor- 
responding to the low eigenvalue behaviour of p(A), we 



obtain the forms 11121) . with a x = (d s - 2)/2 for d s < 4 
while a x — 1 for d s > 4. This expression of a x is analo- 
gous to Eq. J2J), with the Euclidean dimension replaced 
by d s , d^ = 2, djj — 4 and 6 = 1, as appropriate to vec- 
torial models Jg|. This also implies that when d s > d^, 
i.e. T c > |3|, a x > 0, consistently with our previ- 
ous argument. Interestingly, most networks representing 
real physical structures have d s < 2 and, therefore, 
a x < 0. This would imply a diverging response. 

Let us now turn to discrete symmetry systems without 
phase transition, i.e. at d^. In this case T c = and aging 
dynamics can only be observed at T = where Eq. Q 
does not provide a positivity constraint on a x . For in- 
stance, the exact solution of the kinetic Ising chain [J^ 
gives a x = for T — > 0. In the following we intro- 
duce a topological argument showing that the same result 
a x = holds for coarsening systems on FRNs. On these 
structures, each arbitrary large part can be disconnected 
cutting a finite number of links, as those marked x or 
y in the representation of Fig. For FRNs T c = ||, 
and even if the general condition for T c = for discrete 
symmetry models is not known, the available results in- 
dicate that large scale connectivity is the only relevant 
feature and that infinite ramification on scale invariant 
structures implies T c > 0. 

Let us consider phase ordering on a FRN. We assume 
a no bulk flip rule (NBF), where spins in the bulk of 
ordered domains cannot flip. It has been shown 0] 
that this dynamics isolates the aging part of \ leaving 
other properties of the system unchanged. Switching on 
a random magnetic field with expectations hi = 0, hjhj = 
h 2 5ij from time s onwards, Xag can be computed as 

E3 

! N 

X(t,s) = ]mj^22(ai(t)hi), (6) 

i 

where <r, = ±1 and 7TT denotes an average over the 
realizations of hi. Let us focus on a domain of, say, 
up spins whose interface lies at time s on the link de- 
noted x, and is subsequently found in y at time t. The 
domain's size grows from L(s) to L(t). The number 
of spins in the region B between x and i/ is rag oc 
\L{t") — L(s)] d °, where d c is the connectivity dimension 
of the network 2]. We want to estimate the increase 
of the response Xag 9 (A s ) associated to the displace- 
ment of this single interface from x to y. Eq. JHJ) gives 
Xig g {t,s) = -(Nh^-^EhPhfat), where P(y,t) is 
the probability to find the interface in y at time t and 
E h = - J2ieB hi&i = - J2ieB hi is the variation of mag- 
netic energy due to its displacement from x to y. Now we 
make the assumption that the correction to the unper- 
turbed probability Po(y, t) is in the form of a Boltzmann 
factor P h (y, t) = P {y, t) exp(- pE h ) ~ P {y,t)[l - /3E h ] . 
Then x£T ff (M) - -(Nh 2 )-iJ2 y E h [l-j3E h ]P (y,t). 
Since the linear term in Eh vanishes by symmetry we find 
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FIG. 1: Schematic representation of a FRN. 



X#*M = Pm 2 )-^ y E 2 h P (y,t) = 0(Nh*)- l (Ejj). 
Recalling the expression of Eh and the expectations 
of hi one has (E 2 ) ~ h 2 n B . Then P^xtg 9 ^ s ) 00 
A r_1 [L(f) — L(s)] da . The total aging response of the 
system at time t is obtained multiplying Xag 9 {ti s ) by 
the number n,(i) = NL(t)~ dc of interfaces present in 
the system. Then, for t ^> s one has (3~ 1 Xag(t, s) = 
n l {t)(3- 1 Xag g {'t,s) - const. Recalling Eq. ©, this im- 
plies a x — 0. Putting togheter i) and hi) we conclude 
that a x > or a x = is expected in the phase ordering 
of discrete symmetry systems on networks with T c > 
or T c — respectively. 

In order to check this general result, for the Ising model 
we have determined numerically a x in a representative 
class of structures. We have used the NBF to extract 
the aging part of two time functions. Notice that the 
NBF suppresses bulk fluctuations and hence rises the 
critical temperature |l9| . allowing one to avoid numer- 
ically demanding low temperature regions. The para- 
metric plot of Xag(C ag , s) versus C ag , is well suited to 
discriminate between a x = and a x > 0. In fact, since 
C ag (t, s) = git/ s) 0, with g a monotonously decreasing 
function, from Eq. Q one has j) collapse of Xag{C a g, s) 
on a single master-curve, for a x = or jj) lowering of the 
curves increasing s, for a x > 0. 

Coarsening on a percolation cluster is suited to de- 
scribe our results because the bond probability p distin- 
guishes a situation with T c = at the percolation thresh- 
old p = p c [2(|, from those with T c > when p > p c . 
Fig. 12 shows that indeed for the Ising model one has 
a x = at p = p c and a x > for p > p c . Furthermore, we 
have performed analogous simulations on two FRNs with 
T c = 0, the Sierpinski gasket and the T-fractal (Fig. EJ, 
where a x = with great accuracy (Fig. 0J. The situa- 
tion is different for infinite ramification, with T c > 0, e.g . 
on the Sierpinski carpet and the Toblerone lattice |2l| 
a x > (Fig. B} . Similar results are found for different 
discrete symmetry models with non-conserved order pa- 
rameter [(J, such as the 3-state Potts model, confirming 
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FIG. 2: Xag{C a g, s) for the Ising model quenched to T = 1.25 
on the site-percolation cluster at p — p c = 0.407 (upper panel) 
and at p = 0.2 (lower panel). The cluster is built on a 1200 2 
square lattice. Data are averaged over 100 realizations. 



FIG. 4: Xag{C ag , s) for the Ising model at T = 3.0. On FRNs 
a x ~ 0; for the Sierpinski Carpet and the Toblerone lattice 
0.25. 
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FIG. 3: The discrete networks considered in simulations: 
265722 sites Sierpinski gasket, 531442 sites T-fractal, 1048576 
sites Sierpinski carpet and 5078988 sites Toblerone lattice. 

the generality of our result. Interestingly, in all our sim- 
ulations, the value of a x has been found to be extremely 
stable and independent of the temperature at which the 
simulation was performed, even when other asymptotic 
quantities were strongly fluctuating. This suggests that 
a x should be related to a true universal characterization 
of the process. 

In conclusion, we have shown that, in all the cases 
considered, statistical models on networks above dh are 
characterized by a positive non-equilibrium exponent a x , 
while structures at dh have a x — 0. For continuous sym- 
metry systems, our large-TV calculation shows that, be- 
sides this, the whole dependence of a x on the network 



topology can be expressed as for regular lattices, namely 
Eq. , with the spectral dimension d s replacing the Eu- 
clidean dimension d. Our results provide an evidence 
for a relationship between non-equilibrium kinetics and 
large scale topology on general networks and suggest that 
the same topological features of the networks determine 
critical behavior and non-equilibrium exponents of phase 
ordering. 
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